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ABSTRACT 

We present an analysis of the spatial distribution of metal-rich gas in ten galaxy clusters using deep 
observations from the Chandra X-ray Observatory. The brightest cluster galaxies have experienced 
recent AGN activity in the forms of bright radio emission, cavities, and shock fronts embedded in the 
hot atmospheres. The heavy elements are distributed anisotropically and are aligned with the large- 
scale radio and cavity axes. They are apparently being transported from the halo of the brightest 
cluster galaxy into the intracluster medium along large-scale outflows driven by the radio jets. The 
radial ranges of the metal-enriched outflows are found to scale with jet power as i?Fe oc -fjet 42 j with a 
scatter of only 0.5 dex. The heavy elements are transported beyond the extent of the inner cavities 
in all clusters, suggesting this is a long lasting effect sustained over multiple generations of outbursts. 
Black holes in BCGs will likely have difficulty ejecting metal enriched gas beyond 1 Mpc unless their 
masses substantially exceed 10 9 M©. 

Subject headings: galaxies: abundances — galaxies: active — X-rays: galaxies: clusters 



1. INTRODUCTION 

A hot, diffuse plasma composed primarily of hydrogen 
and helium and enriched in heavy elements fills the space 
in and between the galaxies in rich clusters. Chemical 
enrichment of the intracluster medium (ICM) is thought 
to have occurred early, primarily by type II supernova 
explosions. The supernova ejecta raised the heavy ele- 
ment abundance of the ICM to an average of approx- 
imately 1/3 of the solar value (iMushotzkv et al.l 119961 : 
iMushotzkv fc Loewenstemlll997| ). Exceptions are found 
near brightest cluster galaxies (BCGs) located at the cen- 
ters of cool core (cooling flow) clusters. The heavy ele- 
ment abundances there often rise to values approaching 
the mean metallicity of the sun (lAllen fc Fabian! [1998; 
IDupke fc WhTt^mlDe Grandi fc Molendill200lD . Fur- 
thermore, the iron abundance in cool cores is usually 
enhanced with respect to the a-elements. Because the 
a elements are produced primarily by core collapse su- 
pernovae, while iron is produced primarily by type la 
supernovae, the gas in cool cores has probably been en- 
riched relative to the ICM by t he stars associated with 
the brightest clust er galaxies (|De Grandi et al.1 I2004t 
iTamura et al]|200l . 

Several studies have shown that the radial distribu- 
tions of type la ejecta exten d to larger radii than the 
stellar light profiles of BCGs dRebusco et al.| [2005l [20061: 
iDavid fc Nulsenl 120081: iRasera et al.l 120081 ). Assuming 
the ejecta was produced by the stars, the profiles im- 
ply that the ejecta is being displaced from the galax- 
ies by some mechanism. Turbulence induced by cluster 
mergers or outflows fr om active galactic nuclei (AGN) 
are likely possibilities (jSharma et al.ll2009l) . BCGs often 
host radio AGN that are interacting with the hot, metal- 
rich gas by forming cavity s ystems and shock fronts (see 
iMcNamara &: N ulsen 2007, for review). The cavity sys- 
tems in the objects discussed here have displaced between 
10 10 © and 1O 12 of gas. The dynamical forces associated 



with them may also be able to drive the metal-enriched 
gas out of the BCG and into the ICM in outflows ori- 



ented along the rad i o iet (iGopal- Krishna fc Wiita 



2001 



2007 



2003 : lB7uggenl|200a lOmma et al.ll2004t iMoll et al 
iRoediger et al.ll2007l:TPope et al.ll2010D . 

Recent studies of some clusters have identified gas 
with unusually high me tallicity beyond the cor e s and 
along the radio lobes dSimionescu et al.| I2008L 120091 : 
Kirkpatrick et al.1 120095 lO'Sullivan et al.1 120111 ). The 
gas along the radio lobes is, in some cases, enhanced in 
metallicity by up to 0.2 dex relative to the surrounding 
gas. These metal-rich inclusions indicate that the cool, 
metal-rich material originating from the BCG is mixing 
with lower metallicity gas lying at larger radii. It is still 
unclear how general this phenomenon is, and how these 
outflows scale with the properties of the host AGN. A 
systematic analysis of a large sample of X-ray clusters 
with a range in AGN jet power is needed to definitively 
explore this phenomenon. 

Here we show that AGN are driving metals out of 
BCGs to very large distances. We have performed our 
analysis using deep, high resolution X-ray images of 
ten clusters drawn from the Chandra X-ray Observa- 
tory archive. We report on the discovery of a significant 
trend between jet power and the radial extent of up- 
lifted, metal-enriched gas, and we briefly discuss its con- 
sequences. Throughout this letter we assume a ACDM 
cosmology with Ho = 70 km s _1 Mpc -1 , S1m = 0.3, and 
Oa = 0.7. All uncertainties are quoted at the 67% confi- 
dence level. 

2. CLUSTER SAMPLE 

We ex amine 10 cluste r s selec ted fr om the iBirzan et al J 
(|200l . IRaffertv et all (I2006T) and IDiehl et al.1 bops') 
samples of clusters with prominent cavity systems. The 
ten clusters span a broad range of jet (cavity) power. The 
sample includes clusters having Chandra images that are 
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Figure 1. The three panels show the metallicity maps for clusters Abell 262, Hydra A, and Abell 1835, respectively. In all three cases, 
high metallicity gas (shown in yellow) extends outward along the same axis of the cavity system (approximate size and location indicated 
by black circles). The black crosses indicate the approximate cluster centers. 



deep enough to generate metallicity maps with a resolu- 
tion of at least 25 bins with a minimum signal-to-noise 
of 100 per bin. Our goal is to make high precision mea- 
surements of the projected metallicity distribution while 
maintaining a spatial resolution that optimizes the mea- 
surement of the spatial distribution of the high metallic- 
ity gas for comparison to the radio sources. The clus- 
ters that meet these criteria are A133, A262, Perseus, 
MS0735.6+7421, Hydra A, M87, A1835, A2029, A2199, 
and A2597. Each cluster observation was processed us- 
ing version 4.1.2 of CIAO and the calibration database. 
Background flares were excluded using standard filter- 
ing techniques and point sources were identified in the 
combined image of each cluster using wavedetect, and re- 
moved. Blank-sky background files were normalized to 
the source count rate in the 9.5-12 keV band. 

3. ANALYSIS 

3.1. Metallicity Maps and Profiles 

We present example metallicity maps for three clus- 
ters (Abell 262, Hydra A, Abell 1835) in Figure 
1. The combined exposures for each cluster were 
binned using a weighted Voronoi tessellation alg orithm 
(ICappellari fc Copinl [20031: IDiehl fc Statledl2006ri . With 
approximately 10000 net counts or greater per bin, spec- 
tra were extracted and weighted response matrices were 
created using standard CIAO routines. A single temper- 
ature plasma model with absorption (WABSxMEKAL) 
was fit over the energy range 0.5-7 keV. Temperature, 
abundance, and normalization were allowed to vary 
wit h the column dens i ty fro zen at the values quoted 
by IDickev fc Lockm an (1993). The internal ratio be- 
tw een metals were se t to th e solar photospheric values 
of iGrevesse fc Sauvall (|1998l ). Allowing the a-elements 
to vary independently does not show any significant ef- 
fect on the outcome of the iron abundance measurement, 
therefore we have chosen to use the simple model. The 
spectra are reasonably well fit by the single tempera- 
ture model, with the exception of the central region of 
A262. This region is better fit with a multi-temperature 
model which finds a high er central metallicity value (see 
iKirkpatrick et ai]|2009aD . However, this does not affect 
the interpretation of our results as the central regions of 



the clusters are not important for this analysis. 

Metallicity profiles have been measured for each clus- 
ter. Spectra were extracted from regions occupied by 
radio emission and/or cavities in semi- annular shaped 
bins. Each bin size was determined by requiring a mini- 
mum signal-to-nose of 140, which achieves uncertainties 
of approximately 0.05 Zq. A high signal to noise ratio is 
required to find deviations between the profiles along and 
orthogonal to the radio cavity system. These orthogonal 
profiles, which we consider to represent the undisturbed 
region of the cluster, are extracted using the same bin 
width as the radio cavity system profile. The S/N for 
these regions are usually greater than 140 due to larger 
angular sizes. All profiles are fit using the same model 
as the 2D analysis. 

3.2. Pj et -i?Fe Scaling Relation 

The radial profiles for A262, Hydra A, and A1835 are 
presented in Figure 2. The red squares represent the 
metallicity of the radio cavity system and the blue trian- 
gles represent the undisturbed metallicity profile. The 
dotted vertical line indicates the maximum radius at 
which a significant enhancement in metallicity has been 
detected. We refer to this as the iron radius (Rf c ). This 
radius is concisely defined as the location of the radial bin 
furthest from the cluster center where the one sigma er- 
ror bars of the profiles along and orthogonal to the radio 
axis do not overlap. Regions interior to the iron radius 
show a systematic enhancement in iron abundance over 
the undisturbed regions. For regions beyond the iron ra- 
dius, the abundance profiles are indistinguishable. From 
top to bottom the three panels are arranged in order 
of increasing jet power and increasing iron radius. This 
shows that higher power AGN are able to launch metal 
enriched gas to greater distances. 

In Figure 3, the iron radius is plotted against jet 
power derived usin g X-ray cavity data taken from 
iRaffertv et al.l (|2006l ) for all ten clusters in the sample. 
The powers quoted there were derived from the total en- 
ergy of the cavities, ipV, and their buoyant rise time. In 
systems with multiple generations of cavities, the average 
power is used. Note that the iron radius was determined 
using the radial metallicity profiles only and not using 
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Figure 2. Each panel represents metallicity profiles along and 
orthogonal to the radio axis for each cluster. The red square points 
represent the metallicity along the axis of the radio cavity system, 
and the blue triangular points represent the undisturbed region of 
each cluster. The dashed vertical line represents the iron radius. 

the metallicity maps. A trend between jet power and 
iron radius is evident over three decades in jet power. 
The low powered AGN outbursts Pj C t — 10 erg s _1 
drive out material on a scale of a few tens of kiloparsecs. 
For exceptionally large outbursts with jet power exceed- 
ing 10 46 erg s^ 1 , such as MS0735, metals from the core 
are being launched hundreds of kiloparsecs into the ICM. 

To quantify this trend we fit a linear function to the 
logarithms of jet power and iron radius. Performing a 
least squares regression, the best fit plotted in Figure 3 
takes the form 



R Fc = 58 x if e f (kpc), 



(1) 



where jet power is in units of 10 44 ergs s _1 . The RMS 
scatter about the fit is approximately 0.5 dex. The corre- 
lation is strong and shows a fairly small scatter, although 
with only ten objects the true scatter about the mean 
is difficult to evaluate. The scatter will likely increase 
with the inclusion of additional clusters, which will be 
addressed in a future paper. 

A spurious correlation between iron radius and jet 
power, which depends on the volume of the cavities, 
could arise due to the dependence of distance on both the 
linear diameter of the cavity systems and the iron radius. 
This does not appear to be the case. The clusters we con- 
sidered here all have been exposed deeply with Chandra. 
We could have easily detected cool, metal-rich gas on 
all spatial scales of interest here, and at random angles 
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Figure 3. Jet power versus iron radius. The dashed line is the 
best fit to the data. 



with respect to the cavities and radio sources. We did 
not. Hydra A, for example, has a complex cavity system 
that was created by at least three separate outbursts or 
a continuous outb urst that has pers isted for several hun- 
dred million years (|Wise et al.ll2007f ) . The cavity systems 
considered separately have similar jet powers, but their 
centers lie at very different radii between 25 kpc and 250 
kpc. Nevertheless, Hydra A does not deviate significantly 
from the trend in Figure 3, implying that the measured 
iron radius is not a simple function of cavity size alone. In 
Figure 5 we have plotted the inner cavity position for all 
clusters in our sample versus the iron radius. In all cases 
the iron radius lies beyond the inner cavities. This shows 
that the metallicity outflows are not simply tracing the 
current generation of AGN activity, but are maintained 
over multiple generations. The iron radius may provide a 
reliable indicator of average jet power in a cluster where 
cavity power measurements are ambiguous. 

There are uncorrected systematic uncertainties in Fig- 
ure 3. We have made no attempt to measure the addi- 
tional power associated with shock fronts and faint ghost 
cavity systems, which would be difficult to do for the 
entire sample. The total jet power for objects such as 
MS0735 and Hydra A, which have detected shock fronts 
at high significance, are under reported here by roughly a 
factor of two (see McNamara & Nulsen 2007 for a discus- 
sion). Including energy from the shocks effects the slope 
of the fit by ~ 20%. We chose not to boost the power 
of these systems to avoid biasing them with respect to 
other systems. Thus, only cavity power is given here for 
consistency. In addition, the iron radii given here are 
projected values. Their true values may be underesti- 
mated in extreme cases by roughly a factor of two, but 
typically by a few tens of percent. Errors of this size will 
not change our basic result. We intend to explore these 
issues in a future paper. 
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Figure 4. Angle of metal enhanced gas versus cavity angle and 
radio emission angle. The blue squares represent angle on the sky 
with respect to the center of the cavity. Where high resolution data 
is available, the radio angle (330 MHz radio emission) is represented 
by red diamonds. Blue and red points located at the same 6>p e 
highlight the difference in cavity and jet angle in the same cluster. 
The dashed line represents the line of equality. 

3.3. Alignment Between Metal Enhancements and the 
Radio Orientation 

We find that the high metallicity gas outflow regions 
are spatially aligned with the radio and X-ray cavity sys- 
tems. Using the metallicity maps, we have independently 
measured the angle on the sky with respect to the cen- 
ter of the cluster of the high metallicity gas outside of 
the cluster core. We have also measured the mean posi- 
tion angle of the cavity system. The position angle was 
measured using unsharp masks of each cluster in order 
to estimate the locations of the centers of the cavities. 
Finally, the angle of the extended radio emission was de- 
termined for those systems with resolved radio imagery. 
Using 330 MHz images, the radio angle was determined 
by bisecting the radio emission on either side of the nu- 
cleus of the BCG. We have plotted metal angle versus 
cavity center (blue squares) and radio angle (red dia- 
monds) in Figure 4. Both sets of points are consistent 
with the dashed line of equality. In clusters with cal- 
culated radio angles, all three quantities are in excellent 
agreement with one another. 

3.4. Comparison to Independent Metallicity 
Measurements and to Simulations 

Independent reports of metal enriched outflows are 
consistent with the trend shown in Figure 3. For exam- 
ple, the radial abundance effect is seen in a relatively 
low S/N metallicity map of the poor cluster AWM 4 
(jO'Sullivan et all 120101 l20"ll . This study found metal- 
enriched gas extending along the radio jets from the nu- 
cleus of the BCG to a radius of 35 kpc. With a jet power 
of approximately 2 x 10 43 erg s _1 , our scaling relation 
predicts an iron radius of 25 kpc. This result lies within 
one standard deviation of our relationship. 



Figure 5. The radius of the innermost cavity of each cluster is 
plotted against the iron radius. The points represent the cavity 
center and the vertical error bars indicate the radius of the cavity. 
The dotted line of equality is also plotted. 

Recent simulations of AGN outflows have shown that 
radio jets are expected to ad vect ambient mat erial and 
launch it into the ICM (eg., IPope et all 120101). Using 
a 3D AMR hydrodynamical code. iGaspari eTatl (l20ll 
found in their simulations that jets are able to transport 
metal-rich gas from the cores of clusters into the ICM. 
Their simulated jet power of approximately 5.7 x 10 44 erg 
s^ 1 lifted gas 200 kpc from the center before it mixed 
with the surrounding metal-poor gas. The simulation 
over-predicts the iron radius by a factor of 1.6 compared 
to our scaling relation, which represents about a one 
sigma deviation from our relationship. 

Using a series of three dimensional hydrodynamical 
jet sim ulations centere d in a realistic, non-static atmo- 
sphere, iMorsony et al.1 (|2010t ) have evaluated the spa- 
tial impact of radio jets on the hot ICM as a function 
of jet power. They refer to this as the radial "sphere 
of influence" of the jet, which we compare to our iron 
radius. They find that the sphere of influence scales 

1/3 

with jet power as R oc Lj et . This dependence on jet 
power is somewhat shallower than ours. However, the 
rough agreement between our observed scaling relation 
and their simulations is encouraging. 

4. SUMMARY & DISCUSSION 

We have shown that metal-rich gas is being transported 
outward by AGN in all ten clusters in our sample. In each 
case, the metallicity enhancements found in the maps are 
closely aligned with the radio and cavity systems. We 
find a significant trend between jet power and iron radius, 
with a power- law scaling as Rp c oc -Pjc' 42 - This scaling 
relation shows that AGN play an important role in dis- 
persing metals throughout the cores of clusters. A com- 
parison between our scaling relation and recent hydro- 
dynamical simulations discussed in the literature shows 
encouraging agreement. A detailed comparison between 
observed and simulated outflows shows strong potential 
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to place interesting constraints on jet physics. 

Simulations of cluster assembly seem to require out- 
flows in order to spread metals throughout the cores of 
clusters (|Fabjan et al.ll20Tot iBarai et all 120111 ). Major 
mergers at early ti mes are efficient at distributing gas 
throughout the ICM (|Burns et al.ll2008l) , but recent stud- 
ies show in cool core clusters, sloshing due to mergers 
alone cannot explain at later times the large radial ex- 
tent (i.e. > 100 kpc) of mctallicity peaks (Roediger et al. 
2010). A simple extrapolation of our scaling relation 
to a higher jet power indicates that AGN exceeding 
10 47 erg s _1 would be required to drive metals beyond 1 
Mpc, which is a cosmologically interesting regime. Me- 
chanical jet powers of this magnitude have not been 
observed in clusters, although quasars radiating above 
10 47 erg s _1 have been observed in the cen ters of a few 
cool core clusters (eg., iRussell et al.1 [20Toh . Achieving 
mechanical jet powers of this magnitude requires accre- 
tion rates of several tens of solar masses per year onto a 
nuclear black hole. At the same time, mechanically dom- 
inant AGN are expected to occur in radiatively inefficient 
inflow s accreting below a few perc ent of the Eddington 
rate (jNaravan fc McClintocldl2008f ) . Therefore achieving 
such powerful mechanical outflows would be exceedingly 
difficult unless the host galaxies harbor black holes with 
masses significantly in excess of 10 9 M©, which are ex- 
pected to be rare. This may explain why the metallicity 
enhancements at the centers of clusters extend several 
hundred kpc beyond the BCG and generally not to 1 
Mpc. 
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the Natural Science and Engineering Research Council 
of Canada. 
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